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Abstract 

Gravitational waves from inspiralling compact binaries can be reli- 
ably extracted from a noisy detector output only if the template used in 
the detection is a faithful representation of the true signal. In this article 
we suggest a new approach to constructing faithful signal models. 



1 Introduction 

In searching for gravitational waves from an inspiralling compact binary (ICB) 
we are faced with the following data analysis problem: On the one hand, we 
have some exact gravitational wave form h x (t; Xk) where A&, k = 1, . . . ,n\, 
are the parameters of the signal (eg., the masses mi and mi of the members 



retical calculations 
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of the emitting binary) . On the other hand, we ha 
of the motion ofu, and gravitational radiation from 
compact bodies (neutron stars or black holes) giving the post-Newtonian (PN) 
expansions (expansions in powers of v/c) of an energy function E(v) and a 
gravitational wave luminosity function F(v). Here, the dimcnsionless argument 
v is an invariantly defined "velocity" related to the instantaneous gravitational 
wave frequency / GW (= twice the orbital frequency) by v — (jcmf w )3, 
where m = mi + m,2 is the total mass of the binary. Given the energy and 
flux functions one needs to compute the "phasing formula" , i.e. an accurate 
mathematical model for the evolution of the gravitational wave phase <p GW — 
2<f> = F[t; A»], involving the set of parameters {A^} carrying information about 
the emitting binary system. The standard energy-balance equation dE tot /dt = 
—F gives the following parametric representation of the phasing formula: 

where t c and <I> C are integration constants. We now turn to the discussion of 
what is known about the two functions E(v) and F(v) entering the phasing 
formula and how that knowledge can be improved. 
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2 New Energy and Flux Functions 

Let E Tn ee Y,k=o E k{'n)v n and F Tn = Y2=o F k{v)v n , where r) = m 1 m 2 /m 2 
is the symmetric mass ratio, denote the n th -order Taylor approximants of the 
energy -and flux functions. For finite ry, the above Taylor approximants are 
knownB'BO for n < 5. In the test mas&.limit, i.e. 77 — ^ 0, E(v) is known ex. 



the exact flux is known numericallyLru and analytically the flux is known 
up to the order n — 1 1 . The problem is to construct a sequence of approximate 
wave forms h^(t;Xk), starting from the PN expansions of E{v) and F(v). In 
formal terms, any such construction defines a map from the set of the Taylor 
coefficients of E and F into the (functional) space of wave forms. Up to now, 
the literature has only considered the standard map, say T, 

(E Tn ,F Tn )^hl(t,X k ), (2) 

obtained by inserting the successive Taylor approximants into the phasing for- 
mulaOfl. We propose a new map, say "P", based on two essential ingredients: 
(i) the introduction, on theoretical grounds, of two new, supposedly more ba- 
sic and hopefully better behaved, energy-type and flux-type functions, say e(v) 
and f(v), and (ii) the systematic use of Pade approximants (instead of straight- 
forward Taylor expansions) when constructing successive approximants of the 
intermediate functions e(v), f(v). Schematically, our procedure is 113: 

(E Tn ,F Tn ) -> (e Tn J Tn ) -> (ep„,/pj -> (E[e Pn ], F[e Pn , / P J) -» h^(t,X k ). 

(3) 

Our new energy function e(x), where x = v 2 , is constructed out of the 
total relativistic energy E tot using 

e( a )^f ^- m? - m ' ) a -l. (4) 
\ 2mim 2 / 

The function e{x) is symmetric in the two masses. The function E(x) entering 
the phasing formulas is given in terms of e(x) by 

E{x) = [1 + 277 (y/l + e{x) - l)] V2 - 1. (5) 

In the test-mass limit the exact expression for the function e{x) can be 
computed from its definition above which when substitute in Eq. (||) gives the 
well known energy function for a test mass in orbit around a Schwarzschild 
black hole: 

e to = - x T=te< E{x) = -^(i-3^- ^ 
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The test mass exact energy function e(x) has a simple pole singularity while 
the function E(x) has in addition a branch cut. Therefore the function e{x) 
is more suitable in analysing the analytic structure. In the comparable mass 
case, on the grounds of mathematical continuity between the case r\ — > 
and the case of finite rj, one can expect the exact function e(x) to admit a 
simple pole singularity on the real axis cx (x — iTpoic) -1 - We do not know the 
location of this singularity, but Pade approximants are excellent took_for giving 
accurate representations of functions having such pole singularities!!! Indeed, 
it turns out that the Pade_.approximant of the 2PN expansion of e(x) gives 
the exact energy functionE3. We therefore use Pade approximants of e{x) in 
computing the energy function E(x) entering the phasing formula instead of 
the standard Taylor approximants. This greatly improves the accuracy of the 
phasing formula. n 

It has been pointed outB that the flux function F(v; rj = 0) has a simple 
pole at the light ring v 2 — 1/3. The light ring orbit corresponds to a simple pole 
x P oie(v) m the new energy function e(x;rj). Let us define the corresponding 
(invariant) "velocity" i> po ic( r /) = \l x P oic(v)- This motivates the introduction 
of the following "factored" flux function, its Pade approximants fp n , and the 
corresponding flux function entering the phasing formula: 

f(v;r)) = (1- u/wpoic) F(v;t]), F Pri (v;r]) = (1 - v/v po i e )~ fp n (v\rj). (7) 



3 Effectual and Faithful Signal Models 

In order to test whether a given approximant to the wave form is good or not we 
make use of the statistic used in detecting the ICB signal. We shall say that a 
multi-parameter family of approximate wave forms h (t; /J,k), k = 1, . . . , n Mi is 
an effectual model of some exact wave form h x (t; A&); k = 1, . . . , n\ (where one 
allows the number of model parameters to be different from, i.e. in practice, 
strictly smaller than n\) if the overlap, or normalized ambiguity function, 
between h x (t; A&) and the time-translated family h A (t — r; /xj,), 

A(\ k ,Uk) = max — — v y ; v ^" (8) 

r,t ^{h x {t-\),h x {t-\)){h A {t-n)M(t^)) 

is, after maximization on the model parameters /Ltfc, larger than some given 
threshold, e.g. max u , 

A(X k ,fJ,k) > 0.965LJ. [In Eq. (|) the scalar product 
(h, g) denotes the usual Wiener bilinear form involving the noise spectrum 
S n (f)-] While an effectual model may be a precious tool for the successful de- 
tection of a signal, it may do a poor job in estimating the values of the signal 
parameters A^. We shall then say that a family of approximate wave forms 
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h A (t;\£), where the A^ are now supposed to be in correspondence with (at 
least a subset of) the signal parameters, is a faithful model of h x (t; Xk) if the 
ambiguity function A(\k, A;5 ), Eq. (|J), is maximized for values of the model 
parameters A^ which differ from the exact ones Xk only by acceptably small 
biases. A necessary criterion for faithfulness, and one which is very easy to im- 
plement in practice, is that the "diagonal" ambiguity *4(Afc, A;5 = Xk) be larger 
than, say, 0.965. Using this terminology Eq. (||) defines approximants which, 
for practically all values of n we could test, are both more effectual (larger 
overlaps) and more faithful (smaller biases) than the standard approximants 
Eq. (||) . The new sequence of P-approximants exhibit a systematically better 
convergence behavior than the T- approximants t3. The overlaps they achieiic 
at a fixed PN order are usually much higher. From our extensive study Ii3 
of the formal "test-mass limit" r\ = m,\milijn\ + m-i) 2 — > 0, it appears that 
the presently known (u/c) 5 -accurate PN results allow one to construct ap- 
proximants having overlaps larger than 96.5%. Such overlaps are enough to 
guarantee that no more than 10% of signals may remain undetected. 

Table 1: Fraction of events T accessible relative to the case when the true signal is known, 
percentage bias in the estimation of total mass B m = 100(1 — m A /m x ) and percentage bias 
in the estimation of the mass ratio B v = 100(1 — r/ A /r] X ), using wave forms h„ and h^, 
respectively. (A is either T or P, X is for exact, and n is the order of the approximant) 



n 




t p b t b p b t 




Neutron star-black hole binaries 


4 


0.928 


0.998 -6.96 -3.61 12.2 


5.64 


5 


0.945 


1.000 -97.2 -1.11 69.3 


1.83 


6 


0.967 


0.998 1.00 -0.157 -1.56 


0.263 


Black hole-black hole binaries 


4 


0.920 


0.991 1.40 -1.524 0.282 


0.700 


5 


0.532 


0.998 -23.6 -0.205 23.1 


0.042 


6 


0.988 


1.000 0.391 0.019 1.44 


0.119 



Our results are summarized in Table |l|, for two archetypal binaries in- 
vovling neutron stars (NS) and black holes (BH), where we have tabulated 
the fraction of events which the templates constructed out of T- and P- 
approximants would detect relative to the total number of events that would 
have been detectable if we have had access to the true signal. We have also 
listed biases in the measurement of parameters. We clearly notice the superi- 
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ority of the P-approximants. 

Though we believe that the the new approximants are superior over the 
standard ones and shows the practical sufficiency of the presently known 
v 5 -accurate PN results, we still think that it is an important (and challenging) 
task to improve the (finite mass) PN results. Our calculations also suggest that 
knowing E and F to v 6 would further improve the effcctualncss (maximized 
overlap larger than 99.5%) and, more importantly, the faithfulness (diagonal 
overlap larger than 98%) to a level allowing a loss in the number of detectable 
events smaller than 1%, and significantly smaller biases (smaller than 0.5%) 
in the parameter estimations than the present (v/c) 5 results (about 1 — 5%). 
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